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Recent  measurements  of  electron-transfer  (ET)  rates  in  Ru(NHj)5  His  myoglobin  derivatives  have  shown  the 
need  for  improved  theories  for  treating  the  “path"  of  long-range  ET.  We  have  investigated  these  systems 
theoretically  using  a  combined  artificial  intelligence  (Al)-5uperexchange  method.  As  in  our  previous  articles, 
this  model  first  employs  an  AI  search  technique  that  yields  the  important  amino  acid  residues  for  the  media*i(Hi 
of  electrons  between  the  donor  and  the  acceptor.  A  quantum  mechanical  method  is  then  used  to  diagonalize 
the  orbitals  of  the  selected  protein  subset  and  to  calculate  the  electronic  coupling.  Encouraging  agreement  with 
experimental  data  is  found,  no  adjustable  parameters  having  been  introduced.  The  results  yield  the  relevant 
amino  acid  paths. 


iBtrodactkw 


Long-range  electron  transfer  (ET)  is  central  to  many  biological 
processes  such  as  respiration  and  photosynthesis.  Electron- 
transfer  reactions  in  biological  systems  such  as  proteins  have 
therefore  been  the  subject  of  intense  experimental  and  theoretical 
studies.  In  theories  of  electron-transfer  reactions,  factors  which 
influence  the  rate  of  these  reactions  include  the  driving  force 
(-A(7°),  a  nuclear  reorganization  term  (X),  the  separation  distaaae 
(R),  and  the  nature  of  the  medium  separating  the  electron  donor 
and  the  acceptor.*  In  order  to  probe  unambiguously  the  effect 
of  the  intervening  medium  on  the  rate  of  electron- transfer  reactions 
in  proteins,  there  have  been  several  experimental  investigations 
of  intramolecular  electron-transfer  reactions  in  native  and 
modifledproteins.^'’  In  this  paper,  we  present  results  of  theoretical 
studies  on  the  distance  and  medium  dependences  of  the  rates  of 
electron-transfer  reactiems  in  proteins,  particularly  the  ruthenium- 
modified  myoglobin  derivatives  studied  experimentally  fay  Gray 
and  co-workers.' 

For  nonadiabatic  electron-transfer  reactions,  such  as  the  long- 
range  ET  reactions  in  proteins,  the  rate  constant  for  transfer  of 
an  electron  from  a  donor  to  an  acceptor  can  be  expressed  in  terms 
of  a  Golden  Rule  type  expression,  namely,  as  the  product  of  the 
square  of  an  electronic  coupling  matrix  element  (Hoa)  and  a 
nuclear  Franck-Condon  factor  (FC):' 


^Er  =  fl^DAp(FC)  (1) 


The  Franck-Condon  factor  (FC)  can  be  estimated  quantum 
mechanically,  semiclassically,  or  classically.  Various  expressions 
for  FC  are  given,  for  example,  in  a  recent  review.'  In  the  classical 
limit,  the  Franck-Condon  factor  is  given  by 


FC  = 


1 


4XJ?r  / 


(2) 


where  AG”  is  the  standard  free  energy  for  the  electron-transfer 
reaction  at  a  Hxed  donor-acceptor  separation  distance  R  and  X 
is  a  reorganizational  term  which  contains  both  solvent  (protein) 
and  vibrational  contributions. 

The  electronic  matrix  element,  Hda<  describes  the  coupling  of 
the  orbitals  of  the  donor  (D)  with  the  acceptor  (A).  Simple 
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square  barrier  tuimeling  models'^  have  suggested  an  exponential 
decay  of  Hda  with  the  disunce  R  separating  the  donor  and  the 
acceptor: 

"da  =  «°DA  «p{-d("  -  f?®)/2)  (3) 

where  H°oa  the  matrix  element  at  van  der  Waals  contact,  /{** 
(R9  3  A),  and  d  is  the  distance  decay  factor.  Such  an 

exponential  dependence  has  been  experimentally  and  theoretically 
verified  for  several  synthetic  D-A  systems,"-"  where  B  has  been 
found  to  be  between  0.8  and  1.2  A-'.  A  recent  analysis  of  ET 
rates  in  several  biological  systems,  most  particularly  the  bacterial 
photosyntbetic  reaction  centers,  by  Ehttton  and  co-workers  has 
provided  an  overall  support  for  the  exponential  decay  model  for 
Hoa  in  biological  ET  reactions  as  well."  However,  in  electron- 
transfer  systems  which  are  not  structurally  homogenous,  such  as 
jn^oteins,  the  electronic  struaure  of  the  protein  medium  coupling 
D  and  A,  in  addition  to  the  separation  distance,  is  expected  to 
affect  Hoa- 

It  has  been  shown"*  that  the  electronic  coupling  in  derivatives 
of  cytochrome  c  does  not  scale  uniformly  with  distance.  The 
electronic  cou^g  in  these  proteins  is,  cm  the  other  hand, 
satisfactorily  explained  by  simple  pathway  models'*  as  well  as  by 
more  sophisticated  models  which  explicitly  take  into  account  the 
electronic  structure  of  the  protein  medium.''^  The  model  that 
we  have  previously  introduced"^*  to  treat  electr<Huc  interactioDS 
in  proteins  combines  an  artificial  intelligence  (AI)  approach  with 
a  quantum  mechanical  formulation  of  superexchange,  without 
the  introduction  ofany  adjustable  parameters.  The  relative  values 
of  electronic  coupling  elements  obtained  with  this  model  were  in 
good  agreement  with  experimental  results  for  the  cytochrome  c 
derivatives. 

Recently,  Gray  and  co-workers  have  experimentally  studied 
ET  in  Ru-modified  myoglobin  derivatives.'  As  noted  by  these 
authors,  a  simple  pathway  analysis  does  not  adequately  describe 
the  electronic  coupling  in  these  systems.  Nor  did  a  multiple 
pathway  model  which  neglected  interferences  between  pathways. 
Presently,  we  investigate  the  electronic  interactions  in  these 
myoglobin  derivatives  using  our  Al-supcrexchange  method 

Theory 

For  donor-^widge-aoceptor  (D-B-A)  electron-traa^er  systems, 
where  the  distance  between  D  and  A  spans  10-20  A,  the  direct 

®  1993  American  Chemical  Society 


t  Coatrib«itk»  No.  8763. 


6112  The  Journal  of  Physical  Chemistry,  Vol.  97,  No.  23,  1993 


Lettm 


electronic  coupling  between  D  and  A  is  negligible  and  the  protein 
(the  'bridge*}  provides  the  electrtmic  coupling  which  enable  the 
electron-transfer  reaction.  Second-order  perturbation  theory  can 
be  used  to  determine  this  bridge-mediated  electronic  coupling.^-^ 
The  electronic  properties  of  the  bridge,  namely,  the  molecular 
orbitals  and  their  eigenvalues,  are  calculated  separately,  and  the 
interaction  of  the  donor/acceptor  with  the  bridge  is  treated  as 
a  perturbation.  In  a  one-electron  description  of  the  D-B-A 
system,  this  ‘supcrexchange*  mechanism  can  be  formulated  as 


a 


(4) 


where  T.®  (T„^)  represents  the  interaction  of  the  donor  (acceptor) 
orbital  with  the  bridge  orbitals  and  A£a  is  the  energy  of  the 
bridge  orbitals  relative  to  the  energy  of  the  donor  orbital.  The 
above  formulation  includes  bole  transfer  as  well  as  electron 
transfer,  since  the  summation  over  a  runs  over  both  occupied  and 
unoccupied  orbitals  of  the  bridge.  Expressions  for  TJ*.  and 

A£a  can  be  found  in  a  previous  article  (see  eq  8  in  ref  22). 

A  direct  implementation  of  such  a  formulation  for  the  case  of 
long-range  ET  in  proteins  requires  determining  the  electronic 
wave  fimetions  of  the  protein,  which  serves  as  the  bridge  in  such 
reactions.  This  step  necessitates  diagonalizationofa  large  matrix 
(of  the  order  of  6200  X  62(X)  for  myoglobin)  and  is  a  compu¬ 
tationally  formidable  task  and  less  illuminating,  we  believe,  t^  n 
analtmiativeintxedure.  We  have  previously  shown^^  that,  to 
adequatel  j  describe  the  elecmmic  coupling  provided  by  the  protein 
medium,  it  is  not  necessary  to  include  all  the  amino  acid  residues 
of  the  protein  but  only  a  suitably  chosen  subset  of  the  protein. 
This  subset  of  the  protein  is  selected  by  use  of  an  artificial 
intelligencesearch  algorithm.  We  briefly  describe  below  the  way 
the  search  operates  but  refer  the  reader  to  earlier  papcrs“^  for 
a  more  comi^ete  account. 

The  search  starts  from  the  donor  site  O  in  the  protein  and 
attempts  to  locate  the  acceptor  site  A  by  traveling  via  the 
intermediate  aUnns  I  of  the  amino  acid  residues  of  the  protein 
molecule.  These  intermediate  atoms  or  'nodes*  are  chosen  on 
the  basis  of  an  evaluation  function,  EF,  which  is  a  measure  of  the 
coupling  of  the  node  I  to  the  donor  (Fbi)  as  well  as  an  estimate 
of  the  coupling  from  I  to  the  acceptor  (Tia): 


(EF),=  KD,r^A£  (5) 

where  is  a  suitably  chosen  energy  difference.  The  couplings 
Kdi  and  Tu  can  be  estimated  as  explained  in  section  II  of  ref  22. 
It  is  important  to  emphasize  here  that  the  evaluation  function  is 
used  ody  for  the  purposes  of  the  A1  search,  and  the  actual 
electrtmic  coupling  matrix  element  ^da  is  calculated  using  eq  4. 

Using  theevaluati<Hi  function  as  a  heuristic  to  guide  the  search, 
the  AI  search  proceeds  from  the  donor  atom  through  the  protein 
mdecule  to  the  acceptor  atom.  Several  paths,  within  a  preset 
threshold  value  for  the  net  electronic  coupling,  are  determined. 
The  search  is  performed  in  both  directions,  i.e.,  starting  from  the 
donor  and  reaching  the  acceptor,  as  well  as  starting  from  the 
acceptor  and  reaching  the  donor.  The  search  thus  maps  out  a 
number  of  intermediate  protein  atoms  which  enable  it  to  proceed 
frmn  cme  redox  center  to  the  other.  The  amino  acid  residues  to 
which  these  atoms  belong  are  then  considered  as  the  relevant 
amino  acids  for  mediating  the  electron  transfer  in  the  protein. 
The  search  does  not  necessarily  proceed  along  the  backbone  of 
the  protein  but  instead  chooses  the  paths  that  couple  in  the  most 
optimal  way  the  donor  and  the  acceptor  atoms,  ba^  on  both  the 
overlap  of  the  intervening  protein  atomic  orbitals  as  well  as  their 
energies.  Thus,  some  of  the  selected  amino  acids  will  be  bonded 
to  others  not  selected  by  the  search.  These  “dangling*  bonds  are 
completed  by  addition  of  H  atoms.  The  selected  amino  acids 
constitute  the  subset  of  the  protein  that  will  be  used  as  the  bridge 
in  the  calculation  of  Hqa- 


Ru  his  70 


Figure  1.  Amino  acid  residues  determined  by  the  AI  search  for  the  Hit 
70  derivative  of  myoglobin. 


Figwe  2,  Amino  add  residues  determined  by  the  AI  search  for  tbe  His 
48  derivative  of  myoglobin. 

Reaute  aud  Dbeuseinu 

Tbe  proteins  considered  here  are  the  ruthenium-modified 
myoglobin  derivatives  studied  experimentally  by  Gray  and  co- 
workers.  Three  mutant  human  myoglobins  were  modified  by 
replacing  the  heme  by  Zn  mesoporphyrin  (ZnP)  and  coordinating 
a  pentaammineruthenium  complnt  to  a  surface  histidine  residue 
(Ru  His  A'  Mb,  where  =  70, 48,  and  83).  In  the  resulting  singly 
ruthenated  myoglofrin  derivatives,  ET  rates  from  >ZnP*  to  Ru’^ 
were  measured  using  transient  absorption  spectroscopy.*  Tbe 
driving  force  for  tbe  reaction  is  0.82  eV,  and  the  reorganizatiem 
energy,  X,  has  been  estimated  to  be  1.3  eV  from  a  study  of  the 
driving  force  dependence  ET  rates  in  Ru  His  48  Mb.** 
Activationless  ET  rates,  were  extracted  from  the  observed 
ET  rates,  assuming  X  :=  1.3  eV  for  all  tbe  three  derivatives. 

In  tbe  caiculational  procedure,  the  donor  is  the  2b  porphyrin 
group  and  the  acceptor  is  the  ruthenated  histidine  group.  The 
coordinates  of  the  Ru-modified  myoglobins  were  obtained  frmn 
the  refined  protein  oomdinates*'  ^  performing  a  search  of  the 
ruthenated  histidine  conformations.  The  two  side-chain  dihedral 
angles,  C.-C4  and  CfCy,  were  each  rotated  by  1*  fitno  0*  to 
360*  to  find  the  most  stableconformer.**  Tbe  resulting  structure 
files  were  used  in  the  AI  search. 

For  each  of  the  three  derivatives.  His  70  Mb,  His  48  Mb,  and 
His  83  Mb,  the  search  was  allowed  to  select  several  paths,  as 
described  earlier.  A  preset  threshc^d  value  for  the  net  electronic 
coufdingof  any  path  was  used  to  exclude  particularly  unimportant 
pat^.**  Figures  1-3  show  the  amino  acid  residues  that  were 
selected  by  the  search  for  each  of  the  three  myoglobin  derivatives. 
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FIfve  3.  Aouno  acid  residues  determined  by  the  AI  search  for  the  His 
83  derivative  of  myoglobin. 

TABLE  1:  CnlciihiteJ  Ekctnmic  Coopliags  and 
£lcctroii>Trusfcr  Rates 


1  (deriv) 

R(A) 

•ff DA  (cm*')  calc 

calc 

expt* 

1  (His  70) 

11.1 

5.9  X  10-' 

26 

15 

2  (His  48) 

13.8 

2.0  X  10-* 

1 

1 

3  (His  83) 

17.6 

7.0  X  10-< 

0.03 

0.07 

*  Reference  8. 


Only  S-lOof  the  1  S3  possible  amino  acid  residues  are  determined 
by  the  AI  search  as  important  in  these  protein-mediated  ET 
reactions. 

A  quantum  mechanical  calculation  of  the  electronic  coupling 
matrix  element  was  then  performed  using,  as  the  bridge,  the 
amijpo  acid  residues  picked  by  the  AI  search.  As  in  earlier 
papers,^^  extended  Huckel  theory^  was  used  to  calculate  Hda- 
These  electronic  couplings  are  given  in  Table  I.  From  these 
electronic  couplings,  we  can  then  estimate  ik-.-.  the  maximum 
rate  constant  for  ET  in  these  three  derivatives  as  follows;  The 
maximum  rate  constant  for  nonadiabatic  ET  reactions  occurs 
when  the  driving  force  equals  the  reorganization  energy,  i.e., 
when  the  exptmential  factor  in  eq  2  is  unity.  The  maximum  or 
activationless  rate  ccmstant  is  then  given  by  eq  6  (using  a  classical 
ai^noximation): 


The  calculated  values  are  (in  s-')  8.0  X  10^,  1.2  X  10^, 
and  1.1  X  1(F  for  the  His  70  Mb,  His  48  Mb.  and  His  83  Mb 
dmivatives,  respectively.  The  maximum  rates  extracted  from 
the  experimentally  obsmed  rates*  are  7.2  X  10\  3.2  X  10*,  and 
1.8  X 10*  for  these  same  doivatives.  The  good  agreement  between 
these  absolute  values  should  be  regarded  as  fortuitous,  and  it  is 
preferaUe  to  cmisider,  instead,  relatives  values  when  comparing 
calculated  and  miperimental  resultt,  as  in  Table  I;  Results  frcnn 
extended  HOckel  calculations  are  more  reliable  for  relative  values 
of  .^da  *nd  hence  than  for  absolute  values.*'  Further,  the 
‘abstdute*  values  of  km.,  extracted  from  experimental  data  are 
also  dependent  on  the  model  chosen  to  ctHrrect  for  the  nuclear 
factor  and  might  change  somewhat  with  a  different  analysis  of 
the  Rtperimraital  data.  (In  ref  8  and  hm-e,  the  simplest  model, 
nanwly,  the  exponential  expression  in  eq  2,  has  been  chosen,  and 
Exed  X  and  IhG*  have  been  assumed  for  all  the  three  ET  reactions 
in  obtaining  both  the  calculated  and  experimenul  k^  values. 


R  (A) 

Fjgare  4.  Variatiss  of  calculated  and  i—,  wth 

in  the  Ru-modified  myoglobm  derivatives.  The  filled  cudes  refer  to  the 
experimental  values,  the  hollow  circlet  refer  to  the  present  calculatioBS, 
the  filled  triangles  refer  to  the  multiple  pathway  results  of  ref  8,  and  the 
unfilled  triangles  refer  to  their  single  pathway  results.  The  experimental 
and  calculated  km.,  for  the  His  70  Mb  derintive  (R  =  1  i.I  A),  7.2  X 
10*  and  8.0  X  10*  s~',  respectively,  are  too  dose  to  be  distinguisbed  in 
this  plot.  The  k^  obtained  by  t^  single  pathway  analysis  for  His  83 
Mb  (R  =  17.6  A)  is  7.S  x  I0~*  s'*,  which  is  off  this  graph.  As  noted  in 
the  text,  it  is  more  important  to  consider  the  variation  of  these  values 

rather  than  their  absdnte  values. 


The  resulting  “extrapolation*  from  ^et  to  yield  k^  is  not  a  large 
one — a  factor  of  about  S.)  Hence,  it  is  preferable  to  compare  the 
rela  tive  values  of  A—,  for  these  various  derivatives,  and,  therefore, 
in  Table  I  we  give  the  ratios  of  Ab„  with  respect  to  the  His  48 
Mb  derivative.  We  have  also  plotted  the  variation  of  both  the 
calculated  and  experimental  A^u  values  with  distance  in  Figure 
4.  (We  have  chosen  to  use  the  porphyrin  edge-Ru  distances  as 
being  the  most  appropriate  for  the  plot,  since  the  ET  reactions 
soidied  are  *ZnP*  -*  Ru**’  and  in  *2^*  the  extra  electron  resides 
on  the  porphyrin.) 

The  calculated  variation  of  A^u  is  similar  to  the  experimentally 
found  variation,  both  of  them  slwwing  an  approximately  expo¬ 
nential  dependence  with  distance.  However,  the  calculated  values 
exhibit  a  steeper  falloff  with  distance.  In  order  to  test  the 
convergence  of  the  search  procedure,  mtxe  tolerant  cutoff  values 
were  used  in  the  AI  search  enabling  the  selection  of  more  amino 
acid  residues  per  derivative.  The  a^itional  amino  adds  selected 
were  His  64  for  the  His  70  Mb  derivative;  Asp  44,  Leu  49,  Asp 
60,  His  64,  and  His  93  for  the  His  48  Mb  derivative;  and  His  82, 

Ser  92,  and  Pbe  138  for  the  His  83  Mb  derivative.  The  values 

of  Hda  obtained  with  this  more  comprehensive  calculation  were 

6.6  X  10-',  3.6  X  10-*,  and  4.2  X  10^  cm-'  and  the  ratios  of  the 

square  root  of  Aau  were  18, 1,  and  0.01  for  His  70  Mb,  His  48 

Mb,  and  His  83  Mb,  respectively  (normality  with  respea  to  the 

His  48  Mb  derivative).  Comparing  with  the  calculated  ralues 

of  Hda  in  TaUe  I,  it  is  seen  that  there  is  no  dramatic  change  in  ■— 

Hda  due  to  the  addition  of  more  amino  add  residues.  r 

Anuncertaintyinthepresentcalculationistheenagydifreiaice 
in  the  denominator  of  eq  4,  which  depends  on  the  calculated 
value  of  the  energy  of  the  donor  orbital  at  the  transition  state.  j;::* 

In  previous  studies,**-**  we  have  made  use  of  experimental  charge-  U 

transfer  spectra  from  the  donor  spedes  to  the  bridge  to  determine  ^ - — 

this  energy  difference.  Presently,  this  energy  was  set  to  be  the - — 

energy  of  the  LUMO  d'  the  ZnP  group,  since,  as  noted  earlier, 

the  reaction  studied  experimentally  involves  transfer  of  an  electron - - - 

from  the  exdted  sute,  ZnP*.  This  LUMO  is  delocali»d  over  /  ^ 

the  porphyrin  ring,  and,  hence,  extended  HQckel  theory  was  „  Codes 

thought  to  be  adequate  enough  tc  calculate  this  energy'  relative  ^  ^ - 

to  the  energy  of  the  bridge  orWtals,  in  contrast  to  the  case  where or 
the  donor  orbital  is  localized  on  a  metal  center.  This  energy,  lal 
DTIC  QUAUTY  5  j 
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without  further  correction,  was  taken  to  be  the  energy  of  the 
donor  orbital,  which  is  equal  to  that  of  the  acceptor  orbital  in  the 
transition  state.  The  present  work  also  neglects  any  effects  due 
to  conformational  fluctuation  of  the  protein  coordinates  on  the 
electronic  coupling. 

For  comparison,  the  values  of  km  from  the  Beratart-Onuchic 
pathway  analyses  of  these  same  derivatives  from  ref  8  are  also 
plotted  in  Figure  4.  A  B-O  pathway  is  a  single  chain  of  atoms 
connecting  the  donor  and  the  acceptor  and  the  electronic  coupling 
provided  by  a  pathway  is  estimate  by  the  number  of  links  where 
each  bond,  hydrogen  bond  and  through-space  link  contributes  a 
parametrized  decay  factor.  The  two  B-O  pathway  analyses  in 
ref  8  differ,  in  that  one  of  them  uses  a  single  pathway  to  estimate 
the  electronic  coupling  while  the  other  uses  many  such  pathways 
and  sums  the  electronic  couplings  of  the  individual  pathways  to 
yield  the  total  coupling,  neglecting  interferences  between  path¬ 
ways.  Again,  it  is  important  to  note  that  in  the  pathway  analyses 
of  ref  8  only  the  variation  of  km’s  can  be  compared,  since  those 
analyses  do  not  include  the  interaction  between  the  donor 
(acceptor)  with  the  first  (last)  bond  of  the  pathway,  which  can 
cause  a  shift  in  the  absolute  values.  As  was  point^  out  by  the 
authors  of  ref  8,  while  the  inclusion  of  multiple  pathways  improves 
the  agreement,  neither  of  the  two  pathway  analyses  results  in 
satisfactory  agreement  with  the  experimental  data.  Our  procedure 
of  estimating  the  electronic  coupling  automatically  incorporates 
all  possible  amino  acid  paths  within  the  subset  of  protein  selected 
and  further  includes  all  possible  interferences  between  these  paths. 
It  has  been  shown  that  inclusion  of  both  the  interferences  between 
multiple  paths  and  the  non-nearest-neighbor  interactions  are 
important  for  describing  electronic  coupling  in  model  complexes 
and  proteins.^*^^ 

It  is  of  interest  to  compare  the  amino  acid  residues  selected  by 
the  present  AI  search  with  the  (atoms  of  the)  amino  add  residues 
involved  in  the  pioneering  patWay  analyses  of  Beratan  and  co- 
workers.^'  At  present,  only  for  the  His  48  derivative  have  the 
amino  add  residues  from  the  pathway  approach  been  identified. 
For  this  derivative,  the  pathway  analysis  yields  21 1  paths.  All 
of  the  amino  add  residues  involved  in  these  21 1  paths,  except 
one,  are  selected  by  the  AI  search  (see  Figure  2).  Ibe  additional 
amino  add  present  in  the  pathway  treatment  is  Asp  44.  We  find 
that  our  AI  search  proce^  from  Phe  43  directly  to  Phe  46  or 
Lys  45  via  through-space  connections  rather  than  involving  Asp 
44.  It  involves,  only,  in  effect,  2  or  3  amino  add  chain  paths 
instead  of  21 1  individual  atom-by-atom  paths.  When  Asp  44 
was  added  to  the  small  group  of  amino  ad^  depicted  for  the  His 
48  derivative  in  Figure  2,  Hda  changed  only  by  a  small  amount 
from  the  value  of 0.020  cm-*  in  Table  I,  namely  to  0.023  cm.  Asp 
44  did  contribute,  together  with  four  other  amino  adds,  to  the 
larger  set  mentioned  earlier  for  the  His  48  derivative. 

CoBdasioas 

In  the  present  paper,  we  have  studied  the  intramolecular 
electron-transfer  reactions  in  the  ruthenium-modified  myoglobin 
derivatives,  where  recent  rate  measurements  have  shown  the  need 
for  theories  more  sophisticated  than  the  existing  pathway  analyses 
and  homogenous  barrier  models.  We  have  investigated  these 
systems  theoretically  using  a  combined  artifidal  intelligence 
search  and  superexchange  method.  This  model  first  employs  an 
AI  search  technique  that  yields  the  important  amino  adds  for  the 
electron-transfer  reaction  and  then  a  quantum  mechanical  method 
is  used  to  diagonalize  the  orbitals  of  the  selected  protein  subset 
and  to  calculate  the  electronic  matrix  element.  Satisfactory 
agreement  with  experimental  data  is  found.  A  similar  analysis 
has  been  previously  implemented  for  the  ruthenium-modified 
cytochrome  c  proteins  and  was  shown  to  give  results  consistent 
with  experimental  data.^  The  present  approach  can  thus  describe 
the  electronic  coupling  in  both  cytochrome  c  and  myoglobin 
derivatives  without  introdudng  any  adjustable  parameters  and 
offers  a  Hamiltonian-based  method,  which  incorporates  any 


number  of  amino  acid  paths  and  their  potential  interactions  to 
treat  electron  transfer  in  proteins.  We  plan  next  to  investigate 
ET  reaaions  of  particular  biological  relevance  as  in  protein- 
protein  complexes. 
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